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Abstract 

The thermal entanglement in a two-qubit Heisenberg XXZ spin chain is investigated under an 
inhomogeneous magnetic held b. We show that the ground-state entanglement is independent 
of the interaction of ^-component J z . The thermal entanglement at the hxed temperature can be 
enhanced when J z increases. We strictly show that for any temperature T and J z the entanglement 
is symmetric with respect to zero inhomogeneous magnetic field, and the critical inhomogeneous 
magnetic field b c is independent of J z . The critical magnetic field B c increases with the increasing 
\b\ but the maximum entanglement value that the system can arrive becomes smaller. 

PACS numbers: 03.67.Lx; 03.65.Ud; 75.10.Jm; 05.50.+q 



* Corresponding author; electronic address: gfl978zhang2001@yahoo.com 



1 



I. INTRODUCTION 



Entanglement is the most fascinating features of quantum mechanics and plays a central 
role in quantum information processing. In recent years, there has been an ongoing effort 
to characterize qualitatively and quantitatively the entanglement properties of condensed 
matter systems and apply then in quantum communication and information. An important 
emerging field is the quantum entanglement in solid state systems such as spin chains 



. Spin chains are the natural candidates for the realization of the entanglement 



compared with the other physics systems. The spin chains not only have useful applications 
such as the quantum state transfer, but also display rich entanglement features jj. The 



Heisenberg chain, the simplest spin chain, has been used to construct a quantum computer 
and quantum dots . By suitable coding, the Heisenberg interaction alone can be used for 
quantum computation 



12l The thermal entanglement, which differs from the other 



kinds of entanglements by its advantages of stability for the reduction in entanglement of 
an entangled state due to various sources of decoherence and in entanglement in time due 
to thermal interactions are absent as the entanglement at finite temperature takes thermal 
decoherence into account implicitly, requires neither measurement nor controlled switching 
of interactions in the preparing process, and hence becomes an important quantity of systems 



for the purpose of quantum computing. In the studies on t 



re entanglement of Heisenberg 



15 



id ]. It is turned out that the 



spin model, a lot of interesting work have been done w[ 
critical magnetic field B c is increased by introducing the interaction of the ^-component of 
two neighboring spin in Ref [17]. But only the uniform field case is carefully studied in the 
above-mentioned papers. The nonuniform case is rarely taken into account. We know that 
in any solid state construction of qubits, there is always the possibility of inhomogeneous 
zeeman coupling jisl . I19I ] . So it is necessary to consider the entanglement for a nonuniform 
field case. M.Asoudeh and V.Karimipour [2^ studied the effect of inhomogeneous in the 
magnetic field on the thermal entanglement of an isotropic two-qubit XXX spin system. 
We noticed that the entanglement for a XXZ spin model in an nonuniform field has not 
been discussed. Although M.Asoudeh states that the different types of anisotropic 
interactions may not be of much practical relevance to concrete physical realization of qubits, 
in the theoretical analysis we think it is very interesting and should be included in the studies 
of spin chain entanglement. This is the main motivation of this paper. 
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For a system in equilibrium at temperature T, the density matrix is p = 4e - ^ , where 
(3 = -Kp , k is the Boltzman constant and Z = tre~^ is the partition function. For simplic- 
ity, we write k — 1. The entanglement of two qubits can be measured by the concurrence C 
which is written as C = max[0, 2 max M-EjAilU, where A* are the square roots of the 
eigenvalues of the matrix R = pSp*S, p is the density matrix, S = erf <E> erf and * stands for 
the complex conjugate. The concurrence is available, no matter whether p is pure or mixed. 
In case that the state is pure p = with 

= a|0,0) + 6|0,l) + c|l,0) + d|l,l), (1) 

the concurrence is simplified to 

C(if))=2\ad-bc\. (2) 



II. THE MODEL AND THE GROUND-STATE ENTANGLEMENT 

The Hamiltonian of the iV-qubit anisotropic Heisenberg XXZ model under an inhomo- 
geneous magnetic field is 



1 - 



(3) 



i=l 



where J and J z are the real coupling coefficients. The coupling constant J > and J z > 
corresponds to the antiferromagnetic case, J < and J z < the ferromagnetic case. -B > 
is restricted, and the magnetic fields on the two spins have been so parameterized that b 
controls the degree of inhomogeneity. Now, we consider the Hamiltonian for N = 2 case. 
Note that we are working in units so that B, b, J and J z are dimensionless. 

In the standard basis {|1, 1), |1, 0), |0, 1), |0, 0)}, the Hamiltonian can be expressed as 

\ 

J 
J 




/ J Z +2B 
2 



H 



\ 



2 





Jz-2B 



(4) 



A straightforward calculation gives the following eigenstates: 



</>!> = |0,0>, 
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i2> = |M>, 

= ; 1 (-11,0) + 10, 1>), 

1 ,(i|l, ) + |0,l)), (5) 



i + e/j 2 j 



with corresponding energies: 



Ei = ^(J Z -2B), 
E 2 =±(J Z + 2B), 
w Jz 

E 4 =-±+ V . (6) 

where rj = \Jb 2 + J 2 , £ = b — rj and £ = b + r\. Note that when b ^ and J > , the two 
states |03) and [fa) respectively go to the maximally entangled state ^(|0, 1) — [ 1 , 0) ) and 
^(|0, 1) + |1, 0)). For J < , they respectively go to ^=(|0, 1) + |1, 0)) and ^(|0, 1) - |1, 0)). 
We can also find that the eigenenergies are even function of the coupling constant J. So we 
can think the ground-state entanglement exists for both antiferromagnetic and ferromagnetic 
cases and should be symmetric with respect to the coupling constant J. The ground state 
depends on the value of the magnetic field B , the coupling constant J z and r\. It is readily 
found that the ground-state energy is equal to 

E 1 = \{J Z -2B), if V <B- J z - 
E 3 = -lf-r], iirj>B-J z . 

So when rj < B — J z , the ground state is the distangled state \<pi) and when rj > B — J z , 
the ground state is the entangled state \<p3). For each value of the magnetic field B , 
there is a threshold parameter j£ = B — i] above which the ground state will become 
entangled. Accordingly for each value of inhomogeneity rj there is a value of magnetic field 
B-F = i] + J z above which the ground state will loose its entanglement. In the entangled 
phase the entanglement of the ground state is found from (2) and (5) to be 



C(I03» = (8) 



where A = £/J. A = ±1, i.e. 6 = 0, the system enters the maximally entangled phase | ^3) 
with entanglement C(\<t>a)) = 1. This result accords with that in Ref|20|. Here we can also 



know that the ground-state entanglement is independent of the interaction of ^-component 
J z . 



III. THE THERMAL ENTANGLEMENT 



As the thermal fluctuation is introducing into the system, the entangled ground states 
will be mixed with the untangled excited state. This effect will make the entanglement 
decreases. At the same time, the distangled ground state mixes with entangled excited 
states. To see the change of the entanglement, we calculate the entanglement of the thermal 
state p = ' ■ In the standard basis {|1, 1), |1, 0), |0, 1), |0, 0)}, the density matrix of the 

system can be written as 



P12 = ^ 



I e kT o 
e^r(m- 
-s 
\ 



n) 





— s 

e^r [m + n) 




\ 









-Ei 



(9) 



where Z 



bsinh Vr 



k2L . In 



kT (l + ekT) + 2e kT cosh m = cosh n = — 
the following calculation, we will write the Boltzman constant k — 1. From Eq.(9) and the 
definition of concurrence, we can obtain the concurrence at the finite temperature. 

Casel: J z = 0. Our model corresponds to a XX spin model. The eigenvalues and 
eigenvectors can be easily obtained. In Fig.l, we give the results at different temperature 
for the nonuniform magnetic field (B = 0) and the uniform magnetic (6 = 0). From the 
figure, we can know that the entanglement is symmetric with respect to zero magnetic field, 
the nonuniform magnetic field can lead to higher entanglement and double-peak structure. 
This results accord with those seen from Ref 5|. 

Case2: J z = J. In order to compare with the results in Ref [2^, we make the substitution 
J — > 2 J , B — > 2B, b — > 2b. The eigenvalues and eigenvectors can be easily obtained as 



pi) = |0,0) 

¥>2> = |M> 
V3> = 



X 



l + x 2 /J 2 J 



-|1,0> + |0, 1», 



y 



V 1 + v 2 /J 2 J 
5 



mi,o> + |o,i», 



(10) 




FIG. 1: The concurrence for J z = and J = 1 case. T = 0.4 (solid curve) and T = 1.0 (dotted 
curve). The left panel corresponds to the nonuniform case and the right panel corresponds to the 
uniform case. T is plotted in units of the Boltzmann's constant k. And we work in units where B 
and b are dimensionless. 

with corresponding energies: 

e x = J-2B, 

e 2 = J- 2B, 

e 3 = -J -2\J\Vl + 5 2 , 

e 4 = -J + 2\J\Vl + 5 2 . (11) 

where 5 = bj J, x = 1b — 2\J\ + 5 2 and y = 26+Vl + S 2 . For the ferromagnetic case J = — 1, 
the ground-state concurrence is C(\ip3)) = -^rgs • For the ferromagnetic case J = 1, the 
grounc -state concurrence is C(|<£> 4 )) = -yj^t- These results are same with those obtained 
m Ref(20|. In Fig. 2, we give the plot of the thermal concurrence for J z = J = —1 case. 
In order to compare our result with that in Refj3], we let the inhomogeneous magnetic 
field b = 0.458 (this accords to the value of £ in Refj^). We can see that the thermal 
entanglement develops and is maximized for zero magnetic field B and gets the maximum 
value at T = 0. 

CaseS: For any J z . With B = 0, the concurrence as a function of b and T for two value of 
J z are given in Fig. 3. They show that the concurrence are 1 for different J z when 6 = and 
T = 0. At the point, the ground state is |0 3 ) with energy — 4f — 1, which is the maximally 
entangled state and the corresponding concurrences are 1. As the temperature increases, 
the concurrences decreases due to the mixing of other states with the maximally entangled 
state. We can also know that the concurrence decrease with the increasing of From the 
two figures in Fig. 3, we can find that with the increasing J z , the critical temperature T c is 
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FIG. 2: (Color online) The thermal concurrence for J z = J = — 1 case. And the inhomogeneous 
magnetic field b = 0.458. T is plotted in units of the Boltzmann's constant k. And we work in 
units where B and b are dimensionless. 




FIG. 3: (Color online) The concurrence in the XXZ spin model is plotted vs b and T. Coupling 
constant J = 1, the magnetic field B = 0. The left panel corresponds to J z = case and the right 
panel corresponds to J z = 0.9 case. T is plotted in units of the Boltzmann's constant k. And we 
work in units where b is dimensionless. 

improved (for J z = 0, T c is about 2, but for J z = 0.9, T c has a higher value). Thus, we can 
get the higher entanglement at a fixed temperature when J z is increased. 

Fig.4 shows the concurrence at a fixed temperature and magnetic field for three values 
of positive J z . It is shown that the concurrences drop with the increasing value of b and 
arrive to zero at the same b value, which is called the critical inhomogeneous magnetic 
field, for various values of J z . This is to say the critical inhomogeneous magnetic field is 
independent of J z . Moreover, we can see that for a higher value of J z , the system has a 
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FIG. 4: The concurrence in the XXZ model is plotted vs b for various value of J z , where J = 1, 
B = 0.8 and T = 0.6 . From top to bottom, J z equals 0.9, 0.4, 0. T is plotted in units of the 
Boltzmann's constant k. And we work in units where B and b are dimensionless. 




FIG. 5: (Color online)The concurrence in the XXZ spin model is plotted vs T and B, where J = 1, 
J z = 0.4. The left panel corresponds to b = case and the right panel corresponds to b = 0.8 
case. T is plotted in units of the Boltzmann's constant k. And we work in units where B and b 
are dimensionless. 

stronger entanglement which is consistent with Fig. 3. 

In Fig. 5 we give the plot of concurrence as a function of T and B for 6 = and b = 
0.8 respectively when J z = 0.4. For B = and 6 = 0, the maximally entangled state 
|$) = -^(|0, 1) — 1 1 , 0) ) is the ground state with eigenvalue — 4f — \ J\- Then the maximum 
entanglement is at T = 0, i. e. C = 1. As T increases, the concurrence decreases as seen 
from Fig. 5 due to the mixing of other states with the maximally entangled state. For a high 
value of B (in left Fig.5 B = 1.40 and in right Fig. 5 B = 1.70) the state |0i) becomes the 
ground state, which means there is no entanglement at T = 0. However by increasing T, 
the entangled state |03) and |</>4) will mix with the state, which makes the entanglement 
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increase. From the two figures in Fig. 5, we can see that when b is raised, the critical magnetic 
field B c increases, but the maximum entanglement value that the system can arrive becomes 
smaller. 



IV. THE CONCLUSIONS 

In conclusion, we have investigate the thermal entanglement in two qubit Heisenberg 
XXZ spin chain under an inhomogeneous magnetic field. The ground-state entanglement 
and thermal entanglement at a finite temperature are given. We find the entanglement exists 
for both antiferromagnetic and ferromagnetic cases. And the entanglement is enhanced by 
increasing the interaction of ^-component J z . The critical inhomogeneous magnetic field is 
independent of J z . The critical magnetic field B c increases with the increasing |6| but the 
maximum entanglement value that the system can arrive becomes smaller. 
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